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Using samples of 102 million T(IS') and 158 million T(2S) event 
samples collected with the Belle detector, we report on the first search 
for charge-parity-even charmonium and charmonium-like states in T(1S) 
and T(2S) radiative decays. No significant Xcj or r\ c signal is observed 
and 90% C.L. limits on BR(T(nS) -»■ 7x c j) \n = 1,2 and J = 1,2,3) 
are obtained. No significant signal of any charmonium-like state is ob- 
served. The product branching fraction limits BR(T (nS) — > 7X(3872)) 
BR(X(3872) -> tt+tt-J/^), BR(T(n£) ->■ 7 X(3872)) BR(X(3872) ->■ 
n + n-n°J/tP), BR(T(nS) -> 7 X(3915)) BR(X(3915) -> uJ/ip), and BR 
(T(IS') ->■ 7F(4140))BR(F(4140) ->■ </>J/iJ)) (n = 1,2) are obtained at 
the 90% C.L. At the same time, BR(T(25) -> 7X(4350))BR(X(4350) -> 
4>J/ijj)) is also determined at the 90% C.L. Furthermore, no evidence is 
found for excited charmonium states below 4.8 GeV/c 2 . 
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1 Introduction 



Experimental observations near the charm threshold strongly suggest that the spec- 
trum of resonances with hidden charm is remarkably more rich than suggested by the 
standard quark-antiquark template and very likely includes states where the heavy- 
quark cc pair is accompanied by light quarks and/or gluons. Lots of new charmonium- 
like resonances (XYZ particles) in the B factories have been observed in the final 
states with a charmonium and some light hadrons. They could be candidates for 
usual charmonium states, however, there are also lots of strange properties shown 
from these states. These may include exotic states, such as quark-gluon hybrids, me- 
son molecules, and multi-quark states pp. Many of these new states are established in 
a single production mechanism or in a single decay mode only To better understand 
them, it is necessary to search for such states in more production processes and/or 
decay modes. 

States with J PC = 1 can be studied via initial state radiation (ISR) with the 
large T(4S) data samples at BaBar or Belle, or via e + e~ collisions directly at the peak 
energy at, for example, BESIII. For charge-parity-even charmonium states, radiative 
decays of the narrow T states below the open bottom threshold can be examined. 

The production rates of the lowest lying P-wave spin-triplet (Xcj, J=0, 1, or 2) 
and S'-wave spin-singlet (r] c ) states in Y(IS') radiative decays have been calculated 
in Ref. [2], where the former is at the part per million level, and the latter is about 
5 x 10~ 5 . The rates in T(25) decays are estimated to be at the same level [2]. However, 
there are no such calculations or estimations for Ll XYZ particles" due to the limited 
knowledge of their nature. 

The production rates of the P-wave spin-triplet Xcj {J=0, 1> 2) and S- wave spin- 
singlet r] c states in T(IS') radiative decays have been calculated by Gao et al; the 
rates in T(2S) decays are estimated to be at the same level [2]. However, there are no 
such calculations or estimations for U XYZ particles" due to the limited knowledge 
of their nature. 

Belle ever searched for charmonium and charmonium-like states in T(IS') and 
T(2S) radiative decays [3JH]. The data used in those analyses include: (1). a 5.7 fb _1 
data sample collected at the Y(IS') (102 million T(1S) events) and a 1.8 fb _1 data 
sample collected at y/s = 9.43 GeV (continuum data); (2). a 24.7 fb _1 data sample 
collected at the T(2S) peak (158 million T(2S) events) and a 1.7 fb -1 data sample 
collected at ^/s = 9.993 GeV (continum data). 

The numbers of the T(1S) and T(2S) events are determined by counting the 
hadronic events in the data taken at the T peak after subtracting the scaled corre- 
sponding continuum background from the data sample collected at the correspond- 
ing energy point. The selection criteria for hadronic events are validated with the 
off-resonance data by comparing the measured R value (R = ^ e + e ^^+J-)^ ) with 
CLEO's result (5]. 
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In this report, we gave the combined results of search for charmonium and charmonium- 
like states in T(IS') and T(2S) radiative decays in order to compare easily. 

2 Belle and KEKE 

The Belle detector operating at the KEKB asymmetry-energy e + e~ collider [6] is 
described in detail elsewhere [7] . It is a large-solid-angle magnetic spectrometer that 
consists of a silicon vertex detector (SVD), a 50- layer central drift chamber (CDC), 
an array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of 
time-of- flight scintillation counters (TOF), and an electromagnetic calorimeter com- 
prised of CsI(Tl) crystals (ECL) located inside a super- conducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return located outside of the coil is 
instrumented to detect K\ mesons and to identify muons (KLM). 

3 Search for charmonium (-like) states 

We reconstruct J /if) signals from e + e~ or ji + \i~ candidates. In order to reduce the 
effect of bremsstrahlung or final-state radiation, photons detected in the ECL within 
0.05 radians of the original e + or e~ direction are included in the calculation of 
the e + /e~ momentum. In order to improve the J /if) momentum resolution, a mass- 
constrained fit is performed for J /if) signals. Different modes have similar J /if) mass 
resolutions. The J / if> signal region is defined as \M l+i --m J / M ,\ < 30 MeV/c 2 (» 2.5a), 
where mj/^ is the nominal mass of J /if). The J /if) mass sidebands are defined as 
2.959 GeV/c 2 < M t+t - < 3.019 GeV/c 2 and 3.175 GeV/c 2 < M e+e - < 3.235 GeV/c 2 , 
and are twice as wide as the signal region. 

We search for the Xcj in the •yJ/if) mode. The energy deposited by the Xcj photon 
is denoted as 7/ since its energy is much lower than that of 7r (radiative photon 
from T decay). The /i + yU~ mode shows a clear J ftp signal, while the e + e~ mode has 
some residual radiative Bhabha background. Figured] shows the ^iJ /if) invariant mass 
distribution together with the background estimated from the J /if) mass sidebands 
(normalized to the width of the J /if) signal range) for the combined e + e~ and 
modes after all the selection criteria are applied. In Fig. [TJ the left panel is for the 
T(IS') data, where apart from possible weak Xco and Xci signals, the J /if) sideband 
events represent well the signal region, indicating that the production of any of the 
Xcj states is not significant. In Fig. [TJ the right panel is for the T(2S) data, where 
no Xcj signal is observed. There are no structures at higher masses, where we would 
expect excited Xcj states. 

The upper limit on the number (n up ) of signal events at the 90% C.L. is calculated 

f C(x)dx 

by solving the equation = 0.9, where x is the number of signal events, and 

J Q C(x)dx 
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C(x) is the likelihood function depending on x from the fit to the data. 

Bayesian upper limits on the number of events at the 90% CL. are found to be 
11.5, 13.8, and 2.4 for the Xco, Xa, and Xc2 for T(IS') decay, and 2.8, 3.1 and 7.6 for 
the Xco, Xa and Xc2 for T(2S) decay, respectively, 
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Figure 1: The jJ/ip invariant mass distributions in the T(1S) (left panel) and T(2S) 
(right panel) data samples. There is no clear Xcj signal observed. The solid curve is 
the best fit, the dashed curve is the background, and the shaded histogram is from 
the normalized J/ip mass sidebands. 



To study the 77/ c mode, we reconstruct the r\ c mass from the invariant masses of 
K° s K + ti- + c.c, ir+ir-K+K-, 2(K+K~), 2(n + n-), and 3(7r+7r-). Figure [2] shows 
the combined mass distribution for the five r\ c decay modes after all of the selection. 
The left panel is for Y(IS') data, while the right is for T(2S) data. The peak in 
hadronic mass at the J/ip mass, as seen in Fig. |2} can be attributed to the ISR 
process, e + e~ — > JisrJ/iP, while the accumulation of events within the r] c mass region 
is small. The shaded histogram in Fig. [2] is the same distribution for the continuum 
data (not normalized in T(2S) data). 

A simultaneous fit is performed to the five final states. The ratios of the r] c 
(J/ip) yields in all the channels are fixed to BRj6j, where each BRj is the r] c {J/i^>) 
decay branching fraction for the i-th mode reported by the PDG [8], and is the 
MC-determined efficiency for this mode. The fit function contains a BW function con- 
volved with a Gaussian resolution function (its resolution is fixed to 7.9 MeV/c 2 from 
MC simulation) describing the t] c signal shape, another Gaussian function describ- 
ing the J /ip signal shape, and a second-order polynomial describing the background 
shape. The mass and width of the BW function are fixed to the PDG values [5] for 
the rj c . The fitted results are shown in Fig. [21 where the solid line is the sum of 
the best fit functions in the simultaneous fit, and the dashed line is the sum of the 
background functions. 

The fits yield 46 ± 22 and 14 ± 20 T(1S) -> 7^ c and 14 ± 20 T(2S) ^ c signal 
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events, respectively. The corresponding upper limits on the number of the r] c signal 
events are estimated to be 72 and 44, respectively, at the 90% CL. 
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Figure 2: The mass distributions for a sum of the five tj c decay modes from T(1S) 
(left panel) and T(2S) (right panel) data, respectively. The solid curve is a sum of the 
corresponding functions obtained from a simultaneous fit to all the r] c decay modes, 
and the dashed curve is a sum of the background functions from the fit. The shaded 
histogram is a sum of the continuum events (not normalized in T(2S) data). The 
J /if) signal is produced via ISR rather than from a radiative decay of an T resonance. 



The selection criteria for T(1S)/T(2S) -> 7X(3872), X(3872) n+n-J/ip are 
similar to those used for ISR ir + ir~J/if) events in T(4S) data [9]. Except for a few 
residual ISR produced ip(2S) signal events, only a small number of events appear 
above the if)(2S) peak in the tt + tt~ J /if) invariant mass distribution, as shown in Fig. [3] 
(left two plots), where the upper plot is for the Y(IS') data, while the lower is for the 
T(2S) data. For Y(IS'), within the X(3872) signal region, there is only one event 
with a mass of 3.870 GeV/c 2 . The upper limits at 90% CL. on the number of the 
X(3872) signal events are estimated to be 3.9 and 3.6 from T(IS') and T(2S) decays, 
respectively. 

We validate our analysis by measuring the if>(2S) ISR production cross section as 
observed in the 7r + 7r~ J/if> mode. The measured cross section of e + e~ — > •yis-Rip(2S) is 
in agreement with a theoretical calculation of 18.5 pb using PDG jH] values for the 
if)(2S) resonance parameters as input. 

We also search for the X(3872) and X(3915) in the 7r + 7r~7r° J /if) mode. We select 
7r + , 7T~, and J /if) candidates in the X(3872) — > 7i + 7i^J/ip mode and a tt° candidate 
from a pair of photons with invariant mass within 10 MeV/c 2 of the 7r° nominal 
mass. Here the 7r° mass resolution is about 4 MeV/c 2 from MC simulation. Figure [3] 
(right two plots) shows the 7i + 7i^7i°J/ip invariant mass distributions, where the open 
histogram is the MC expectation for the X(3872) signal plotted with an arbitrary 
normalization. We observe two events in the n + tt~ it J / if) mass spectrum between 
3.6 GeV/c 2 and 4.8 GeV/c 2 in the T(1S) data. For these two events, the 7r + 7r - 7r°J/^ 
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masses are 3.67 GeV/c 2 and 4.23 GeV/c 2 , and the corresponding 7r + 7r tt° masses are 
0.54 GeV/c 2 and 1.04 GeV/c 2 , respectively. The event at 3.67 GeV/c 2 , is likely to be 
from e + e~ — > 7isr^<//'0 — > 7isR 7r+7r ~7r°£ + £~, since 0.9 events are expected from MC 
simulation. In the T(25) data, there are a few events scattering in the 7r + 7r~7r°J/ , 
mass spectrum. No event is observed within the X(3872) or X(3915) mass region 
in the T(1<S) data, while there is one event with m{ii + n~ -k q J / ip) at 3.923 GeV/c 2 
and m(7r + 7r~7r °) at 0.790 GeV/c 2 from the T(2S) data. We determine n up for the 
number of X(3872) signal events to be 2.3 and 4.2 at the 90% C.L. for T(15) and 
T(2S) decays, respectively. Assuming that the number of background events is zero, 
the upper limits n up for the number of X(3915) signal events are 2.3 and 4.4 at the 
90% C.L. for T(IS') and T(2S) decays, respectively. 

We search for the K(4140) and the X(4350) (only in the T(2S) data) in the (frJ/ip 
mode. The selection criteria are very similar to those in the analysis of X(3872) — > 
7i + 7i^J/ip and the 4> is reconstructed from a K + K~ pair. According to MC simulation, 
the signal region is defined as 1.01 GeV/c 2 < M K+K - < 1.03 GeV/c 2 . After 
applying all of the above event selection criteria, there is no clear J/ip or <fi signal 
in T(IS') or T(2S) data samples. Nor are there candidate events in the F(4140) or 
X(4350) mass regions. The upper limits on the number of V(4140) and X(4350) 
signal events are both 2.3 at the 90% C.L. 

4 Summary 

To summarize, we find no significant signals for the Xcj or r] c , as well as for the 
X(3872), A(3915), F(4140), or X(4350) in T(1S) and T(2S) radiative decays. In 
addition, we find no evidence for excited charmonium states in the invariant mass 
distributions of all final states below 4.8 GeV/c 2 . Table [1] lists the final results for 
the upper limits on the branching fractions of all the states studied. In order to 
calculate conservative upper limits on these branching fractions, the efficiencies have 
been lowered by a factor of 1 — cr sys , here a sys is the total systematic error. The 
results obtained on the Xcj and ^ c production rates are not in contradiction with the 
calculations in Ref. 0. No A(3872), X(3915), or F(4140) signals are observed, and 
the production rates of the -n + n~J/ip, n + n~ n J / tp , uJ/ip, or 4>J/ip modes are found 
to be less than a few times 10 -6 at the 90% C.L. 
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Figure 3: (a) Distributions of the 7r + 7r~ J /ip invariant mass, (b) Distributions of the 
7r + 7r~7r° J/ip invariant mass. The upper two plots are for T(IS') decay, while the 
lower for T(25) decay. Points with error bars are data, open histograms are the MC 
expectation for the X(3872) signal (not normalized). The peak at 3.686 GeV/c 2 in 
(a) is due to i[>(2S) production via ISR. 



Table 1: Summary of the limits on T(IS') and T(2S) radiative decays to charmonium 
and charmonium-like states R. Here BR(T — > 7_R) up (£>#) is the upper limit at the 
90% C.L. on the decay branching fraction in the charmonium state case, and on the 
product branching fraction in the case of a charmonium-like state. 
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